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Abstract 



We present an exact solution of the anisotropics of cosmic background radiation 
(CBR) from a local collapse described by a spherical over-dense region embedded in 
a flat universe, with the emphasis on the relationship between the dipole (AT/T)^ 
and the quadrupole (AT/T)^ anisotropy. This result has been used to examine the 
kinematic quadrupole correction (AT/T)^ = (AT/T)^/2, which is usually applied 
to remove the contamination of the quadrupole produced by local density inhomo- 
geneities when finding the cosmic amplitude of the quadrupole at the surface of last 
scattering. We have found that the quadrupole of local collapse origin cannot al- 
ways be approximately described by the kinematic quadrupole. Our numerical result 
shows that the difference between the kinematic and local quadrupoles depends on 
the size and matter density in the peculiar field, and the position of the observer. For 
a given dipole, the local quadrupole can be different from the kinematic quadrupole 
by a factor as large as 3. Therefore, the kinematic quadrupole correction remains 
an uncertain factor in the determination of the amplitude of a cosmic quadrupole. 
Nevertheless, a preliminary analysis shows that this uncertainty might not dramati- 
cally change the cosmological origin of the COBE-DMR's quadrupole, unless a huge 
pecuhar gravitational field is assumed. 

Subject headings: cosmic background radiation - cosmology: theory 
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1. INTRODUCTION 



This paper is aimed at studying the contribution of local collapse of matter to 
the anisotropies of cosmic background radiation (CBR), especially the relationship 
between the dipole and the kinematic quadrupole. 

It is generally believed that the CBR dipole anisotropy comes from a Doppler 
effect. If the observer moves with velocity v with respect to the CBR rest frame, the 
special relativity (SR) Doppler effect would lead to a frequency-independent thermo- 
dynamic temperature distribution on CBR (Peebles & Wilkinson, 1968) 



T = To 



V 1 s T 

1 + - COS * + -(-) COS 2* + 
c 2 c 



(1) 



The terms in the right-hand side are, respectively, the monopole, the dipole and the 
kinematic quadrupole anisotropies. 

Eq.(l) is often used 1) to determine the peculiar velocity v, and 2) to calculate 
the correction of the kinematic quadrupole. For instance, the recent COBE-DMR first 
year sky maps show a dipole amplitude (AT)^ = 3.365 ± 0.027 mK toward direction 
= (264.4° ±0.3°, 48.4° ±0.5°) (Smoot et al, 1992; Bennett et al, 1992; Kogut 
et al, 1993). If the entirely observed dipole results from our peculiar motion, the 
kinematic quadrupole anisotropy should be (AT)^ = (AT)^/2T ~ 2.1/iK. This result 
has been used by the COBE-DMR team to obtain the CBR quadrupole anisotropy 
(Smoot et al, 1991, 1992). They removed the kinematic quadrupole term from the 
DMR maps in order to eliminate the influence of local density inhomogeneities on 
the quadrupole amplitude. Since the amplitude of the kinematic quadrupole is about 
13% of the cosmic quadrupole, the difference between the quadrupoles produced by 
the SR kinematics and local density inhomogeneities would be one of the factors 
leading to uncertainty in the amplitude of the cosmological quadrupole. Therefore, it 
is necessary to study the conditions, under which the kinematic quadrupole correction 
given by eq.(l) is valid. 

The initial velocity of cosmic matter underwent a decrease due to the expansion 
of the universe. It can then be reasonably assumed that our present peculiar velocity 
results totally from the infall motion toward the center of the local collapse of matter, 
i.e., the peculiar motion of the observer is completely given by the gravitation of 
local matter clustering. In linear approximation, the relationship between the matter 
density fluctuation S and the pecuhar velocity is (Peebles, 1980) 

v^^{HoD)n^-^5 (2) 

where D is the distance to the center of the perturbation, Q, the density parameter 
and Hq, the Hubble constant. Thus, eqs.(l) and (2) imply that the CBR dipole 
and kinematic quadrupole anisotropies are essentially caused by the local density 
inhomogeneities. In other word, the CBR's dipole and quadruple anisotropy should 
be calculated as an effect of a locally time-dependent density inhomogeneity. 

The question is then: can the contribution of a local collapse of matter to the 
CBR dipole and quadrupole anisotropies be approximately described by a Doppler 



3 



effect of eq.(l), i.e. a pure SR effect? Obviously, the SR Doppler effect description 
is principally different from that of a local collapse. The former is a kinematic effect, 
of which the only parameter is the observer's velocity v. The latter is a dynamical 
model, and governed by the matter distribution or the initial density contrast 5, the 
location of the observer and the size of the local collapse. 

Eqs.(l) and (2) also tell us that the dipole anisotropy depends on the first order of 
the density fluctuation 5 and then the kinematic quadrupole is of the second order. It 
is well known from general relativity that in a linear approximation the behavior of a 
comoving object in an expansion or collapsing metric can be equivalently described as 
a Doppler motion. But such an equivalence will no longer be held if the higher orders 
are involved. Indeed, in terms of the dipole term one can not distinguish between a SR 
Doppler effect and a local collapse. The velocity v determined by a dipole anisotropy 
is just equal to the v caused by the locally gravitational collapse. However, one can 
not expect that the SR kinematic quadrupole remains the same as the quadrupole 
caused by the local collapse. 

The original intention of the kinematic quadrupole correction in the COBE sky 
maps is to remove the contamination of a quadrupole given by the local density 
inhomogeneity. Therefore, in principle, this correction should not be done by the 
kinematic quadrupole, but the local collapse of matter. We need then to conduct a 
quantitative comparison of the kinematic quadrupole with that of a local collapse. 

The CBR anisotropies from a collapse at redshift greater than 1 have been ad- 
dressed by several authors (e.g. Rees & Sciama, 1968; Dyer, 1976; Olson & Silk; 
1979; Raine & Thomas, 1981; Kaiser, 1982; Occhionero, Santangelo & Vittorio, 1983; 
Nottale, 1984; Dyer & Ip, 1988; Arnau et al, 1993). However, the contribution from a 
local collapse to the CBR anisotropies, especially the relationship between the dipole 
and local quadrupole, has not been carefully studied thus far. We will investigate the 
CBR dipole and quadrupole anisotropies from a locally spherical collapse described 
by a Tolman-Bondi universe. This model has recently been used to analyze the CBR 
anisotropy produced by perturbations with sizes of 100 ~ 1000 Mpc (Fang and Wu, 
1993, hereafter Paper I). The advantage of this model is that one can find all needed 
exact solutions, which allow us to check the availability of the kinematic quadrupole 
correction of eq.(l). 

In Section 2, we will discuss the metric and null geodesic in a spherically sym- 
metric collapse. Section 3 shows the linear, the high-order and the exact numerical 
solutions of the dipole and quadrupole anisotropies. In Section 4, a relationship 
between peculiar velocity and CBR anisotropies will be presented. Finally, a brief 
summary and discussion will be given in Section 5. 

2. METRIC AND NULL GEODESIC 

For simplicity, we consider a spherical density perturbation embedded in a flat 
universe. The universe can then be generally modeled as a Tolman-Bondi metric 
(Paper I): 

ds'^ = e^(^'*)(ia;2 + r\x, t){de^ + sin^ edcj)^) - df (3) 
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and 



1 + xmm - ^ ^ 

' ^ Pci ' 

where H^, p{x,ti) and pd are, respectively, the Hubble constant, the density and the 
critical density of the universe at epoch U. We denote that ' — d/dx and d/dt. 

The dynamics of the collapse from the initial perturbation at x = is described 
by the collapse factor S{x, t) which is defined as 

r^S{x,t)x, S{x,ti)^l (5) 

The dynamical equation of S{x, t) is 

Pci 'J Pci 

Let the initial density perturbation be 5{x), the density distribution p{x, ti) in eqs.(4) 
and (6) can then be written as 

■p{x,ti)= pci{l + S{x)) (7) 

If the initial density perturbation 5q is assumed to be constant in the range of x < x^ 
one has 

§(rj.\ ^ J X <Xc /gV 

\ ^o{Xc/xf, X> Xc. 

In this case, the dynamical equation (6) can be solved analytically (Paper I). 
The Oth-component of the null geodesic in the metric of eq.(3) is 

where is the Oth-componcnt of a photon's four momentum and a is an affine 
parameter of the null geodesic. From the condition of ds = 0, one can find the energy 
shift of a photon, which is assumed to be emitted at t = te with frequency u^. and 
received at {xo,to) with frequency uq, 

The trajectories of the photon can be obtained from 



dt \ \dt 



(11) 
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Considering the trajectory of a photon is approaching a straight hne when r is large, 
one can obtain d(l)/dt from solving the eq.(12), and then the energy shift eq.(lO) and 
the trajectory eq.(ll) become 



£ = e/"" '^^'^"^'^ {l - j^" UUTj (13) 

^ _ ±g-A/2 ( ^ C \ 

dT \ l-J^'^UtedTj ^ ^ 

where the new time and space coordinates (T, X) arc defined by T = t/tc and 
X = x/te- Therefore, Tq = to/te, Xq — Xo/tg and Xc — Xc/tg. U and ^ in eqs.(13) 
and (14) are, respectively, 

U-n^-"-) (15) 

^=(^)%in^*e^^^*^'^^ (16) 

with Tq = r{xo, to), and \1/ is the incidence angle of the photon, i.e. the angle between 
the directions of the photon's trajectory and the line of sight from observer to the 
center of the perturbation. 

Since the universe is flat, the CBR anisotropy produced by the local collapse is 
given by 



AT 



where Tq = (1 + ZdY^"^ and za is the redshift at decoupling time (te). Without a 
loss of generality, we will take ti — te, i.e., the S{x) is the density fluctuation at the 
recombination time t^- 

Because we are interested in the effect of a local collapse, in the following, only 
the case of Xq < x^. or Xq < X^ will be considered. Using the expansion of S{X, T) 
given in the Appendix, one find the solution of AT/T up to the second order of Sq as 

^ = - ASdT + XI sin^ * r' -^,5dT 



Tco „ AX 



SdT 
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3/1 (Tco) + Xq sm ^ 



+2dov\Q sm W 



rpV3 ;^(0)2 ^2/3 ^(0)2 



+ 2X0' sin^ * / — / /i5dr + X^ sin^ * / ^^5ciT 



(18) 



where and -Uj (i = 1, 2) are defined in the Appendix. AX is the first-order correction 
to the photon trajectory, and ATc, the cross time correction when the photon enters 
into the perturbation regime. They are, respectively. 



AX = ± 



(X(o))2 _ sin^ ^ rTo F(X(°), T) 
Jt 



X(o) 



2^2/3 



-dT 



JTrXI 



F(X(o),T) 



(19) 



(20) 



where X^^^ is the zero-order solution of the photon's trajectory, which can be found 
from eq.(14) by taking 5 = 0. It is 

(X(o))2 = Xl sin^ * + (STo'^^ - 3rV3 _ ^^2 (2I) 

Similarly, the zero-order solution to the photon cross time (Tc) is found to be 



T. 



cO 



(22) 



The function F{X,T) in eqs.(19) and (20) is defined as 



X2 - X^ sin^ * 
S^{To)^ S,{T) rT 



g^^o - + /^^ A^dT + Xl sin^ ^ |L,,T 



(23) 



and Qi is also given in the Appendix. 



3. SOLUTIONS OF AT/T 
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3.1 First- Order Solution 



The first-order solution of AT/T can be obtained by substituting the zero-order 
solution of the photon trajectory X^^^ of eq.(21) into the first two terms of eq.(18). 
After a straightforward computation, the first-order solution is found to be 

.| + At.'/»X.c„s*-4^.0(T„-)) (24) 

The above expression is also an expansion with respect to the parameter (I/Tq). The 
largest term is of the order of (I/Tq)^^/^. Eq.(24) shows that in the approximation 
up to the first-order of So, AT/T consists mainly of two parts: a monopole term 




'AT\ fX^ Xi 2 X3 



T \15 45 135 To 



So (25) 



and a dipole term 



(l^) -^^o^'^o^ocos* (26) 

The physical explanations of these results are simple. The monopole is given 
by the gravitational redshift of the local matter inhomogeneity So, which leads to an 
isotropic increase of the CBR temperature. The dipole term depends on the distance 
between the observer and the center of the density perturbation (Xq). In the case of 
Xo — 0, i.e., the observer sits at the center of the local collapse, the dipole term will 
disappear. This indicates that the dipole anisotropy is indeed due to the asymmetry 
of the local collapse around the observer. However, it should be pointed out that 
the dipole anisotropy depends on the time-dependence of the local inhomogeneity. If 
the gravitation field around the observer is static, the asymmetry [Xq ^ 0) do not 
produce such a dipole anisotropy. 

Comparing eq.(26) with eq.(l), we have 

V = ^T^^'XoSo (27) 

This means that the CBR dipole anisotropy produced by a local collapse can be 
equivalently described as a SR Doppler effect if the observer is assumed to have a 
peculiar velocity given by eq.(27). Let's show that eq.(27) is indeed the observer's 
peculiar velocity produced by the gravitation field of the local collapse. The proper 
distance D corresponding to X — x/te is 

D^- — ^^X (28) 
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where Zi is the redshift when the perturbation occurred. We will take it to be the 
redshift at decoupling era, i.e. Zi = z^. In the paper I, we have found that for a given 
^0 the present density contrast of the local collapse should be 

^^\{l + za)5^. (29) 

In a flat universe, p — Pcr — SHq/SttG. Using eqs.(28) and (29), eq.(27) can be 
rewritten as 

„ = |(«„^)(^)._^, (30) 

where g = GAM/D^ is the observer's acceleration raised by the extra-mass AM = 
(47r/3)L>^Apo- Therefore, v in eq.(27) is the same as that in eq.(2). In a word, in the 
linear approximation it is reasonable to determine the observer's peculiar velocity by 
the CBR dipolc anisotropy and the Doppler effect formula (1). 

If the CBR dipolc is totally given by the local collapse, one can then find from 

eq.(26) that Xo6o = {15/2){AT/T)dTo^^^ ~ 2.6 x 10"^. This is, in fact, a constrain to 
the local collapse causing the dipole. For instance, if we assume the size of this local 
collapse is of the order of the distance to the Great Aractor, i.e. Xq ~ 1, the initial 
density perturbation 6q in this area should be about 5o ~ 2 x 10~^, or from eq.(29) 
today's density contrast is about 1 x 10~^. 

3.2 Second- Order Solutions and Quadrupole Anisotropy 

In a similar way, the second-order solution of AT/T can be found from eq.(18) 
as follows 



AT 



(31) 



which is also written in the scries of (1/To) up to the order of (l/To)-2/3. Eq.(31) 
indicates that, up to the 6q approximation, the dominant term is still the dipole, 
because the amplitude of dipole is of the order of Tq, while the amplitudes of the 

monopole and quadrupole terms are only of the order of Tq^^. 
The quadrupole anisotropy of the local collapse now is 

{^) = i-^T^^'^yoCOs2^ (32) 



Comparing this amplitude with that in eq.(26), one find 

(33) 



'AT\ ^ 



T 




9 



This is just the SR relationship between the anisotropies of the dipole and kinematic 
quadrupole. Therefore, the SR kinematic quadrupole correction is available till the 

approximation of eq.(32) is correct. However, when the terms of the order of Tq^^'^, Tq 
are taken into account, the dipole-quadrupole relation of eq(33) will no longer hold 

true. For instance, up to the order of 5q and Tq^^, eq.(33) should be corrected as 



Tn 



1/3 



(34) 



where the factor is given by 



19^c 
3780 



1 



140X 



+ 



229^3 
61440X3 



+ 



81920X5 



+ 



4096XJ 



41Xn 



1333X3 



+ 



467X0^ 



3833Xo^ 



9800X, 2064384X3 5734400X5 



+ " +0(^1 

11468800XJ ^Xr 



(35) 



Eq. (35) shows that the term of the SR kinematic quadrupole does not always dominate 
the quadrupole anisotropy shown in eq.(34). Therefore, the SR kinematic quadrupole 
may not always be a good approximation of the quadrupole produced by a local 
collapse. In principle, the kinematic quadrupole correction of eq.(33) should also be 
replaced by the local quadrupole correction of eq.(34). 

Let's define a ratio between the local quadrupole to the SR kinematic quadrupole 



as 



(AT/T), 



(1/2)(AT/T)i 



(36) 



Figure 1 plots the dependence of q on X^ and Xq as given by eq.(35). One can see 
that q does not sensitively depend on the observer distance but the whole size of the 
local inhomogeneity. The difference between local and kinematic quadrupole, i.e. that 
q significantly deviates from 1, mainly occurs in two cases: 1) very small clustering 
with a scale less than a few Mpc, and 2) very large clustering with a scale greater 
than 103 Mpc. 



3.3 Numerical Solution 



In order to accurately compare the kinematic quadrupole with the local collapse 
model, we made the numerical solutions of the dipole and quadrupole amplitude of 
the local collapse model. These solutions depend on three parameters: 1) the initial 
density fluctuation do, 2) the size of local collapse Xc, and 3) the distance of the 
observer to the center of the local collapse Xq. 

We still lack the detailed knowledge of the local collapse. It is generally beheved 
that the bulk motion of horizon sized volume is negligible. Therefore, the Local 
Group's peculiar velocity, which may govern the dipole term, was probably induced 
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by the local inhomogeneities in the density field with size comparable with horizon at 
the time of decoupling. We have shown in Paper I that the COBE-DMR result of CBR 
anisotropy on an angular scale of 10° implies the existence of collapses on scale as large 
as about 1000 h^Q Mpc with an initial density fluctuation 5o ~ (2.8 — 6.9) x 10~^, or 
its present density enhancement is (1.7 — 4.2) x 10~^. On the other hand, the distance 
of our Local Group to the center of the collapse should at least be greater than the 
distance to the Great Attractor, which is estimated to be 80 hgo^Mpc (Lynden-Bell et 
al, 1988), i.e. Xq — 0.7. Therefore, it would be valuable to consider the following two 
cases: Xc = 1.4 (~ 150 hso^Mpc) and 10 (~ 1000 h^g'^Mpc), i.e., the lower value of Xc 
is about 2 times of the distance to the Great Attractor and the higher value of Xc is 
about the size of horizon. 

The numerical results are listed in Table 1, in which [exact] indicates the exact 
solution of the dipole amplitude, and [Sq] and [S^] denote, respectively, the solutions 
up to the first- and second-order of Sq. The density fiuctuations at the decoupling 
epoch (zd = 10^) are taken to be 10~^, 10^^ and 10~^, respectively. The corresponding 
values of the SR kinematic quadrupole (1/2)(AT/T)^ = (1/2) [exact] ^ have also been 
listed for comparisons. The term of (AT/T)^ is the exact solution of the quadrupole 
amplitude. 

One can find from Table 1 that: a) For the dipole anisotropy the linear approx- 
imation is good if the initial density fluctuation is less than lO^'^. b) The dipole is 
nearly independent of the size of the local collapse, c) The local quadrupole sensitively 
depends on the size of the perturbation, d) The kinematic quadrupole is always larger 
than the exact solution of the local quadrupole. e) The difference between the kine- 
matic and exact quadrupoles does not vanish with the decrease of density perturbation 
Sq. In the case of Xc — 10 the ratio 1/q can be as large as 3. 

4. DIPOLE AS THE FUNCTIONS OF DENSITY CONTRAST AND 

PECULIAR VELOCITY 

In this section, we intend to represent the CBR's dipole anisotropy by the ob- 
servable parameters such as the present density contrast and peculiar velocity. Up to 
the second order of 5o the dipole anisotropy is (see eqs.(24) and (31)) 



First, the present density contrast Ap/p as a function of initial fluctuation 6o is 
given by 



Using eq.(Al-A5), one can find the expression of Ap/p expanded as a series of Sq 




(37) 




(38) 




3^2/3 ,2 
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2625" 



175^° 



34 



^1/3 

-in 



92 



-1/3 



375 ° 175 ° 125 V 



875 ° 525" 



+ 21 ° 



Reversing this expansion and retaining the first three terms, we have 

2 



1 (^P 



5 85 / A// 
3 " 63 I V 



14075 (/\p 
11907 V V 



(39) 



(40) 



Substituting this result into eq.(37), one find the expression of dipole up to the second- 
order of the present density contract Ap/p, 



(41) 



Therefore, in hnear approximation, the dipole anisotropy would be overestimated from 
the measurement of the local density fluctuation. 

Second, the peculiar velocity inside the perturbation regime can be derived from 
the continuity equation (Peebles, 1980) 



dAp/p 
dT 



So 



V ■ (1 + Ap/p)v = 0. 



(42) 



It is 



1 ,^ ./3 f dAp/p \ So 
Using eqs.(38), (39) and (40), the pecuhar velocity can be expanded by Ap/p as 



(43) 



(HoD) 



1 (Ap; 



Ap' 



4 (Ap\ 328 
63 \ p~) + 11907 y~p 



+ ■ 



(44) 



This result is consistent with the recent work by Gramann (1993), who showed that the 
peculiar velocity would be overestimated by the linear approximation. Substituting 
this expression into eq.(37), we have the solution of the dipole anisotropy as a function 
of peculiar velocity 



'AT^ 



. 25 / c 
1 + — ' 



63 \DHo 



(45) 
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or 



V /AT\ , / 25 Ap' 



(46) 



Therefore, the SR relation eq.(l) can be used to determine pecuUar velocity from the 
CBR measurement only if the local density contrast Ap/p is very small. Considering 
the fact that the present density contrast Ap/p is of the order of 1 on the scale of 
superclusters, the higher order correction for the peculiar field may not be negligible. 
The peculiar velocity divergence in our model is simply 



or 



V- V 



21 V P y 3969 \ p 1 



(47) 



(48) 



1 + 0.190(Ap/p) - 0.046(Ap/p)2 + ... 
A similar result has been empirically found recently by Nusser and Dekel (1991, 1993). 



5. DISCUSSION AND CONCLUSIONS 



It is usually believed that in linear approximation the CBR anisotropics consist 
primarily of three parts, namely, 1) a Doppler effect from the motion of the observer 
with respect to the CBR rest frame, 2) a Sachs- Wolfe effect at the surface of last 

scattering and 3) a time-dependent potential effect along the photon path (Martinez- 
Gonzalez, Sanz & Silk, 1990). In this description, one can not, in fact, distinguish 
between the effects of a timc-dcpcndcnt potential and a Doppler motion. The contri- 
bution of a local collapse of matter to CBR anisotropy can totally be contained in the 
Doppler effect of the observer's peculiar velocity. One can then use the linear rela- 
tionship between the peculiar velocity v and local density contrast Ap/p to calculate 
the CBR dipole anisotropy (AT/T)^ [cq.(l)]. 

However, in terms of the CBR quadrupole anisotropy, the equivalence between 
the Doppler effect and the local gravitational field will no longer exists. This is be- 
cause kinematic quadrupole is essentially non- linear. As a result, the quadrupole 
anisotropy produced by a local inhomogeneity is not generally equal to the SR kine- 
matic quadrupole. Especially, if the peculiar field is comparable with the horizon size, 
the relationship between the dipole and the quadrupole terms will be substantially 
different from that given by the SR Doppler effect. We have found that, based on 
a simply spherical density perturbation model, the local quadrupole can be different 
from the kinematic quadrupole by a factor as large as 3. 

The purpose of the SR kinematic quadrupole correction used in the data reduction 
of COBE observation is to remove the local quadrupole amplitude from the CBR sky 
temperature maps, and then the remaining maps will contain only the component of 
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the cosmic quadrupole, i.e. the component given by the density fluctuation on the 
last scattering surface. The kinematic quadrupole correction is, in fact, the correction 
of local quadrupole. Yet, as we have shown, the local quadrupole can not be simply 
replaced by the SR kinematic quadrupole. Therefore, the original purpose of the 
kinematic quadrupole correction may not be completely achieved by removing the 
quadrupole given by eq.(l). 

The quadrupole amplitude of the local density inhomogeneities depends on the 
matter distribution, the density contrast, the size of the local gravitational field and 
the position of the observer. All these parameters are poorly known at the moment. 
Therefore, it seems to be very hard to precisely identify the local quadrupole. A 
possible way to obtain the local quadrupole would be to analyze the temperature 
map to see if they contain a component with polar axis paralleled to the direction 
of dipole. This is because the axes of both the dipole and local quadrupole should 
point out to the center of the local inhomogeneity. Principally, without a precise 
local quadrupole correction, the observed quadrupole should be regarded as a sum of 
cosmological and local components. Thus, the local quadrupole correction leads to 
an uncertain factor in finding the cosmological quadrupole term from observed CBR 
maps (e.g., the COBE measurement). 

In spite of the fact that the local quadrupole may be remarkably different from 
the kinematic quadrupole, our simple model shows that for a given dipole amplitude 
the amplitude of the local quadrupole is always less than the corresponding kinematic 
quadrupole, if the local peculiar field is of the order of 100 Mpc in size. This means 
that the kinematic quadrupole seems to be an upper limit to the local quadrupole. 
Because the presently observed quadrupole amplitude by COBE is about 8 times 
larger than this upper limit, one may concludes that the uncertainty of the local 
quadrupole would not lead to an uncertainty greater than 13% of the observed cosmic 
quadrupole. 

Certainly, this conclusion is model-dependent. The present model is actually a 
toy one, which has assumed the simplest matter distribution - a constant density 
contrast. It may largely deviate from the real local matter distribution. For instance, 
the peculiar velocity distribution derived from this toy model [eq.(44)] shows a linear 
increase with the distance from the perturbation center, which does not fit with the 
observations around the Great Attractor (Faber & Burstein, 1989). Moreover, if the 
initial density perturbation is assumed to be highly non-spherical, the local quadrupole 
may have an amphtude greater than (AT/T)^. Therefore, the loophole hidden in 
the kinematic quadrupole correction has not been totally closed yet. One needs to 
simulate this effect by using a more realistic model of local collapse of matter, e.g., 
the models given by N-body simulations. The study of the interaction between the 
local matter distribution and the CBR quadrupole would be of great significance for a 
better understanding of both the initial perturbation and the local collapse of matter. 

The authors thanks Dr. J. Heine for his helps. WXP is grateful to the CNRS 
and the Wong K.C. Foundation for financial support. This work was also partially 
supported by NSF contract INT 9301805. 
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APPENDIX 

For the initial perturbation eq.(8), the collapse factor S{x,t) can be expanded in 
the series of Sq as 

SiX,T) = So{T) + S,{T)5o + S2{T)5l + S,iT)6l + ... (Al) 
where the coefficients Si {i = 0, 1, 2, • • •) can be determined by eq.(6) to be 

So{T) = 7"/=^ (A2) 

o T^2/3 7^4/3 

^■(^) = -T5rvs + — -— 

1 4 42-1/3 2^2/3 2^4/3 37.2 

S^iT) - -2257^4/3 + 315^1/3 + 9~ ^ T5~ ~ 175 ^ ' 



SsiT) 



11 2 2 



10125TV3 4725rV3 1125r2/3 945rV3 1575 

5J.2/3 2TV3 23r«/=^ , 

H ^ \ (A5) 

81 875 45 175 7875 ^ ^ 

Similarly, one can find the expansions of r, A, ^ and U with respect to 5o by their 
definitions (Paper I) 

2r' 

A = TT m 

and 

as well as eqs.(15) and (16). Let 

^ = ^ + /i'^o + « + ... (A8) 

^ - - = U160 + ^2^0 + - (A9) 
z r 

e"'^' = ^ + ^1^0 + ^2^0 + - (AlO) 
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where the coefficients fi — fi{T), Ui — Ui{T) and gi — gi{X,T) {i — 1,2,3, •• •) are 
found to be 
1. for X < Xc, 

2 2 

2_2_2 4_ 26T1/3 

/2 - ~ 352-2 ~ 225T4/3 ^ 45TV3 525 ^ 

xi- = 0, i = 1,2,3,- •• (A13) 

_ 1 2 1 2^2 

~ 5 15rV3 ~ 3r2/3 ~ 9T2/3 *^ 

_ 1 1 32 2 2T2/3 



5 45r8/3 315T5/3 9T2/3 35 
_ 2X2 4^2 ^ 2X^ 

~ ^ ~ 135r5/3 + 27T2/3 ~ 81T2/3 ^ ^ 

2. for X > Xc, 

/i - + (A16) 

^ ^8 2_ 4TV3 

■^^ ~ 7573 35T2 + 225T4/3 45TV3 +525 ^ ^ 

-1 = - + ^ (A18) 

_ 2 2 2 2T1/3 

«2 - + - + (A19) 

_ 2 4 ^ 2X2 

" ~5 " 15T57^ ^ 3T27^ ~ 92^ ^ ^ 

1 11 92 12 1 



5 225^8/3 315^5/3 25r 9^2/3 

13T2/3 4X2 gj^2 4J^2 2X 



4 



175 45 135rV3 27^2/3 81^2/3 



(A21) 
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Table 1 Comparison of the SR Kinematic Quadrupole with the Exact Solutions 



size 


So 


(AT/T), 


[(AT/T),]V2 


(AT/T), 


[So]/So 




[exact] /(5o 


Xc = 1.4 


10-4 
10-5 


2.951 
2.951 
2.951 


3.457 
3.002 
2.957 


3.638 
3.003 
2.957 


6.6 
4.5 
4.4 


6.0 
4.1 
4.1 




10-3 


2.951 


3.457 


3.632 


6.6 


2.6 


Xe= 10 


10-4 


2.951 


3.002 


3.000 


4.5 


1.4 




10-5 


2.951 


2.957 


2.953 


4.4 


1.4 
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Figure Captions 



Figure 1. The ratio q of the local quadrupole (AT/T)g to the SR kinematic quadrupole 

(1/2)(AT/T)^. Only the terms with orders greater than T^^^ in eqs.(34) are consid- 
ered. This plot illustrates the dependence of the ratio q on the size of the inhomo- 
geneity and the position of the observer. For the collapse with intermediate scale of 
about 10^ Mpc, we have g a; 1. For both the very small and very large scale collapses, 
the local quadrupole will significantly be different from the kinematic quadrupole. 
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